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We study numerically a Ginzburg-Landau type equation for micelles in two dimensions. The 
domain size and the interface length of a cellular structure are controlled by two feedback terms. The 
deformation and the successive splitting of the cellular structure are observed when the controlled 
0^ ' interface length is increased. The splitting instability is further investigated using coupled mode 

' equations to understand the bifurcation structure. 
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Complicated chemical reaction in a confined cellular structure is considered to be an important step to life. Oparin 
considered that " Coacervate" played an important role in the origin of a cell in prebiotic chemical evolution [l| . The 
' compartmentation by some membrane structure is important for a cell to be independent of the external atmosphere. 
We studied the creation and reproduction of model cells with semi-permeable membrane 0. Vesicles and micelles 
^ can take a cellular form and are considered to be a model system of primitive cells 0] . Self-replication of reverse 

I \ micelles Q and vesicles [1, by the increase of the number of the membrane molecules were observed in experiments. 
O Various types of deformation of the cellular structure such as budding, splitting and birthing were observed. 

On the other hand, the control of spatio-temporal patterns has been an important topic of nonlinear dynamics. 
The spiral patterns and the spatio-temporal chaos were controlled by some feedback mechanisms [1, Q . We studied 
' the control of domain size in the Ginzburg-Landau type equation and the method was applied to the problem of cell 
I differentiati on fiol . In this paper, we try to control the interface length in a Ginzburg-Landau type model for 
micelles [l2l.ll3l| in two dimensions. We will find a splitting instability of cellular domains in the model. Although the 
, splitting instability of two-dimensional puls es was observed in the numerical simulation of the Gray-Scott model (T^ 
^ ■ and an experiment of the FIS reaction [l5| | , it is important from a viewpoint of artificial cells to study the splitting 
instability of the cellular structure in such a micelle model. 
Our analysis is based on a free energy functional: 
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^ ■ i^[0(r)] = / dr[{c/2){V^f + (l/2)g(V0)2 - (1/2)02 ^ (^/4)^4 _ ^ {I / l2)b{V cj>)% (1) 

in : 

' where </> denotes an order-parameter such as the difference of oil and water concentrations in a problem of the mixture 
of oil, water and surfactant. The region with a large value of |V(/)| corresponds to the interface region including 

\^ ' the surfactant. The surface energy is controlled by parameters g and b. The larger area (length) of the interface is 

' preferable in case of negative large values of g and b. A time-dependent Ginzburg-Landau equation is given by 

^' ^ = -^ = 0-</,3^^_ cV^^ + g^^^ + 6(V</.)2V20. (2) 

> We consider the control of 5i = / drcj) and 5*2 = / drlVcj))'^ in two dimensions. Ifc = /i = & = and g > 0, there is a 

5h ' domain wall solution 



4> = ta,nh{x/y'2g). (3) 

Then, is proportional to the domain-size difference of domains satisfying (j) — 1 and domains satisfying (j) — —1, 
and S2 is approximately expressed as S2 ~ 4:/{3y/2g)l, where I is the total length of the interface between the two 
domains. The control of Si and S2 to certain fixed values leads to the control of the domain size and the interface 
length. If c,/i and b are not zero, the above approximation is not always good, but we call the control of S2 the control 
of the interface length in this paper. We control 5*1 and ^2 to certain fixed values by changing the parameters /i and 
g using the global negative feedback as 

^ 7('S'io - -Si), ^ = 7('S'2 - S'20), (4) 

where 7 is a decay constant for fi and g, and Sio and ^20 are target values of 5*1 and 5*2. If Si is larger (smaller) than 
the target value 6*10, fi decreases (increases), which leads to decrease (increase) Si. Similarly, 5*2 is larger (smaller) 
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FIG. 1: Cellular domains satisfying </> > at (a) 520 = 50, (b) 55, (c) 65, and (d) 70. 
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FIG. 2: Cellular domains satisfying </> > at (a) = 75, (b) 100 and (c) 275. 



than the target value 5*20, g increases (decreases) and the interface region decreases (increases), which leads to decrease 
(increase) 5*2. As a result of the negative feedback effect, Si and ^2 are expected to approach 6*10 and 5*20. 

We have performed numerical simulation using the pseudo spectral method with Fourier modes of 128 x 128. The 
system size is L x L = 40 x 40. We have assumed the integration range for Si and S2 as a circular region of radius 
L/2 = 20 for the comparison with the analysis of coupled mode equations below. We have obtained qualitatively 
similar results even if the integration range is assumed to be the total square region of L x L. The target value 6*20 is 
increased slowly from S'20 = 25 to 325. That is, 6*20 is stepwise increased as S20 = 25 + n at t = lOOn where n is an 
integer. The initial condition is assumed to be </) = 1 inside of the shghtly elliptic region (x— L/2)^ + 1.05(?/ — L/2)^ = 9, 
and (j> = —1 outside of the region. Figures l(a),(b),(c) and (d) show cellular domains, where > is satisfied, at 
S20 = 50,55,65 and 70. The domain takes a circular form at ^20 = 50. The deformation to an elliptic form begins 
at S20 ~ 53, which leads to a dumbbell shape at 5*20 — 65 and finally it is split into two domains at S20 ~ 68. The 
split pattern is observed at 520 — 70 in Fig. 1(d). When 520 is further increased, more cellular domains appear by 
the deformation and the splitting. Figure 2(a),(b) and (c) shows the cellular domains respectively at 6*20 — 100, 175 
and 275. The number of cells increases stepwise as 3,4 and 5. 

The deformation and the splitting instability of a circular domain can be approximately expressed by coupled mode 
equations. From the direct numerical simulation of Eqs. (2) and (4), it is expected that 0(r) is approximated at 
4>o{r, t) + (j)2{r, t) cos(26'), where r and 9 are polar coordinates around the center {L/2, L/2). The substitution of the 
approximation into Eqs. (2) and (4) yields 



dt 



dt 




4^ , (5) 
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FIG. 3: (a) Profile of (t>{x,y) (solid curve) at y = L/2 for = -1500 and = 50 in Eqs. (2) and (4). Profile of (l>{x) by 
Eq. (8) (r is set to be x.) and the mirror image for a; < (dashed curve) at the same parameters. The two curves are well 
overlapped, (b) Eigenvalue for Eq. (9) as a function of 520 for 5io = —1500. (c) Mean amplitude S3 of the perturbation (ji'iir) 
as a function of ^20 at 5io = —1500 for Eqs. (5), (6) and (7). 



dr dr \ dr'^ r dr J \ dr"^ r dr J \ dr^ r dr^ dr^ dr J 



and 



dg 

dt 




L/2 

dr2Trr(j)a{r) 



dr{(dcj,^ldrY + {l/2){d^2/drY + - ^20 ■ (7) 



Note that 02 must behave 02 ~ a^r'^ + a^r^ + ■ ■ ■ near r — because the angle dependence is cos(26'). There exists 
always a solution with the circular symmetry satisfying 02 = 0. Such a circular solution 00(7") obeys the equations: 



— = 0^ f (3/2)0„0v I II I hj^^^X ( I ^'^"^ \ g( I 
dt \ dr J \ dr"^ r dr J \ dr'^ r dr 

_ j^SVo , 2d'<j>o 1 9^00 , 1 900 

H 

d9_ 

dt 



7 I S'lo - / dr27rr0o(r) I , 











7 / dr(90o/ar)^ - ^20 . (8) 



These equations are obtained from Eqs. (5) and (7) by setting 02 to be zero. The linear stability of the circular 
solution can be investigated by the linear equation obtained from Eq. (6): 



'2 _ , Q i,2 , , I, "(Po 



9^02 , 1 502 402 \ , ^,d4>od4)2 f d^4>o , 1900 



dt 't'o't'^ + \ dr ) \ dr^ ~^ dr dr r^ ) ^ dr dr \ dr"^ ' r dr 

'd^(t>2 1902 402\ fd^4>2 2 5^02 992^2 9 902 



(9) 



The solid curve in Fig. 3(a) shows the profile of (f>{x^y) at y — L/2 at Sio — —1500 and 5*20 — 50 for Eqs. (2) 
and (4). The dashed curve in Fig. 3(a) is a stationary solution to Eq. (8). The two curves are well overlapped, and 
the difference is hardly seen. That is, the approximation by Eq. (8) is good. Figure 3(b) shows the eigenvalue of 
the linear Eq. (9) as a function of ^20 for a fixed value of 5io ~ —1500. The instability occurs at 6*20 — 52, which 
is consistent with the critical value S20 ~ 53 by the direct numerical simulation of Eqs. (2) and (4). Figure 3(c) 
displays 5*3 = {Jq^"^ drr(j)2}^^^ as a function of 6*20 obtained numerically for Eqs. (5), (6) and (7). It means that 
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FIG. 4: (a) (/f>o(0) as a function of for 5*10 = —1500. (b)Profile of (j)o{x) (solid curve) and the mirror image for x < 
at 5*20 = 85 (solid curve) and 5*20 ~ 80 (dashed curve), (c) Deformed domain at ^20 = 80. In the shaded domain, 
00 (r) + (p2{r) cos(20) > 0. (d) Split domains at 520 = 85. In the shaded two domains, fpoir) + </'2(r) cos(20) > 0. 



the supercritical bifurcation occurs at S20 ~ 52. That is, the eUiptic deformation grows continuously. The splitting 
instability is also approximately described by Eqs. (5), (6) and (7). Figure 4(a) displays the value 00(0) at r = as a 
function of S20 for Sio = —1500 obtained by numerical simulation of Eqs. (5), (6) and (7). A discontinuous transition 
occurs at 820 ^ 82. The profile of <j>o{r) has a peak at r = for 820 < 82. On the other hand, (j>o{r) has a peak at 
nonzero r for 820 > 82 as shovirn in Fig. 4(b). Figure 4(b) displays the profiles of (poix) at 6*20 — 80 and 85. The 
discontinuous transition of the profile 00 (^) is clearly seen. Figures 4(c) and (d) show the deformation of the cellular 
domain at (c) S20 = 80 and (d) S20 — 85. In the shaded regions, 0o(^) + 02(^)cos(20) > 0. The two-peak structure 
shown in Fig. 4(b) appears as a two-cell structure in Fig. 4(d). The critical value 5*20 ~ 82 in the coupled mode 
equations is larger than the critical value S'20 ~ 68 of the splitting instability by the direct numerical simulation by 
Eqs. (2) and (4). It is partly because the higher modes including cos(2m6') with to > 2 is truncated in Eqs. (5), (6) 
and (7). 

In summary, we have proposed a Ginzburg-Landau type model for micelles under the control of the domain size and 
the interface length. As the interface length is increased, a circular cell is deformed to an elliptic form and then split 
into two cells. By increasing further the interface length, many cells are created by the deformation and the splitting 
instability. We have proposed coupled two-mode equations and found that there are two successive bifurcations for the 
splitting instability. One is the supercritical bifurcation, where the circular symmetry is broken continuously. At the 
second bifurcation point, the splitting of the cellular structure occurs discontinuously. In the problem of micelles, we 
can interpret that the increase of the interface length corresponds to the increase of the surfactant materials created 
by some chemical reactions inside of the micelles. The splitting processes might be interpreted to correspond to the 
self- replication process of micelles found in the experiments 0]. 
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